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Abstract 

This paper considers N mobile nodes that move together in the vicinity of each other, whose initial poses as 
well as subsequent movements must be accurately tracked in real time with the assist of A/(> 3) reference nodes. 
By engaging the neighboring mobile nodes in a simple but effective cooperation, and by exploiting both the time- 
of-arrival (TOA) information (between mobile nodes and reference nodes) and the received-signal-strength (RSS) 
information (between mobile nodes), an effective new localization strategy, termed "cooperative TOA and RSS" 
(COTAR), is developed. An optimal maximum likelihood detector is first formulated, followed by the derivation 
of a low-complexity iterative approach that can practically achieve the Cramer-Rao lower bound. Instead of using 
' simplified channel models as in many previous studies, a sophisticated and reaUstic channel model is used, which can 

(N : effectively account for the critical fact that the direct path is not necessarily the strongest path. Extensive simulations 

bX)! are conducted in static and mobile settings, and various practical issues and system parameters are evaluated. It is 

^ ■ shown that COTAR significantly outperforms the existing strategies, achieving a localization accuracy of only a few 

\ tenths of a meter in clear environments and a couple of meters in heavily obstructed environments. 

Index Terms 

. mobile localization, wireless localization, user cooperation, maximum Ukelihood estimation, iterative detection, 

' time-of-arrival, received-signal-strength 
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I. Introduction 

Efficient and accurate localization is desirable for a variety of wireless ad-hoc networks, especially mobile ad 
^ hoc networks and robot systems that will be forming dynamic coalitions for important missions like search-and- 
CN rescue and disaster recovery [1]-[14]. While the world-wise deployment of global positioning systems (GPS) has 
^ significantly alleviated the localization problem [3], there exist several scenarios, such as inside a building or in a 
^ forest, where GPS signal is lacking. Further, because of the hardware cost, size, battery and other concerns, not all 

the nodes in a system may be equipped with a GPS receiver. 
O Wireless localization is being studied in the system architecture, protocol, and algorithm level. This paper 
O considers mobile localization in the algorithm level, and attacks position traclcing, the fundamental computation 
. . problem at the heart of mobile localization. In some situations, the initial pose of the mobile node is known and 
. ^ hence only incremental errors need to be compensated. The more challenging case, also referred to as the "global 
^ localization problem," is when the initial position and all of the subsequent movements must be detected from the 
scratch [5]. A desirable mobile localization algorithm should provide good accuracy and robustness. The former is 
particularly relevant to indoor applications where localization errors should in general be controlled below a few of 
meters or even a few tenths of a meter. The latter is pivotal to the system survivability in a dynamic environment 
where the communication channels may be randomly faded, shadowed, or heavily obstructed without line-of-sight 
(LOS). Finally, the good quality of localization should be achieved in a cost-effective manner with simple hardware, 
since many mobile nodes, especially those in large systems, may be constrained by their cost and sizes, and hence 
are only provisioned with limited resources like a single small omnidirectional antenna, an inexpensive low-precision 
clock and a small piece of bandwidth. This paper aims to develop new localization strategies that can achieve all 
of these goals. 

A fundamental method to determine the point location of a target node (i.e. a robot or a cell phone) is triangulation. 
Existing techniques exploit different modalities of the radio frequency (RF) signal, such as received angle-of-arrival 
(AOA), received signal strength (RSS), time-of-flight (TOE) (also known as time-of-arrival or TOA), and time- 
difference-of-arrival (TDOA), to compute and estimate the location or the relative location of the target node with 
respect to the reference nodes or anclior nodes[l], [11]. 
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OCT1133027. Bao and Li are with the department of electrical and computer engineering, Lehigh University, Bethlehem, PA 18015. Email: 
{xib3, jingli}@ece. lehigh.edu. 



Since the AOA technique requires the availability of expensive directional antennas at the reference nodes and/or 
at the sensor nodes, this is much less employed in practical systems than the other techniques. The RSS technique 
exploits the received signal power together with a path-loss and shadowing model to provide a distance estimate 
between the sender and the receiver. It works effectively when the nodes ai^e close to each other, but is quite sensitive 
to shadowing, multi-path fading and scattering, and hence the performance deteriorates quickly with the increase 
of the distance [7]. The TOA technique uses the signal propagation time between a target node to a reference 
node to compute the target position in a circular track. The TDOA technique estimates the distance difference 
of a target node to two reference nodes through received-time difference, thus determining the position of the 
target in an ellipse track without knowledge of the signal transmitting time. Comparing to the RSS technique, 
the TOA/TDOA-based methods can achieve a reasonable accuracy at a much larger measurement range. It is also 
possible to improve the accuracy of TOA/TDOA by increasing the signal transmit power and the system bandwidth, 
but since the localization inaccuracy is largely caused by multi-path and no-line-of-sight (NLOS), the TOA/TDOA 
accuracy will hit a bottleneck at some point, beyond which further increasing transmit power or bandwidth does 
not help [14]. To combat NLOS and other channel dynamics and to improve positioning accuracy, hybrid schemes 
that employ multiple different modalities are proposed [9] [10]. For estimating the incremental pose change of a 
mobile node, a posteriori methods, including Kalman filters, multi-hypothesis Kaknan filters, Markov localization 
and Monte Carlo location [2] [5], are also shown to be useful. 

This paper considers the global localization problem in a highly dynamic environment where the communication 
channels may be randomly faded, shadowed, or heavily obstructed without line-of-sight (LOS). One big challenge 
about timing an RF signal in an obstructed environment is that the direct path signal is not necessarily the strongest 
path, and hence the detector tends to miss the direct path signal but catch the strongest path signal, causing the 
measurement of the propagation time to be longer than it actually i^ However, the majority of the existing studies 
have only looked into clear-sight environments, and hence have not explicitly accounted for this critical source of 
error. In this paper, we will adopt a sophisticated and realistic 2-ray channel model which is derived from real-world 
experiments and which can effectively capture the reflective path distortion [14], in addition to all the other usual 
noise and measurement error. 

Unlike the majority of the work that considers localization as independent tasks of individual (mobile) nodes, 
here we consider cooperative localization. Many applications involve locating and tracking multiple target nodes 
that stay in the vicinity of each other. Examples include tracking first-responders operating in a buddy system in a 
search-and-rescue task in a building, positioning a group of scientists exploring a forest, or monitoring a robot that 
is equipped with two (or more) sensors on both sides of the body performing tasks in a cave. Exploring practical 
application needs, we propose to engage neighboring nodes in a simple but effective cooperation, and to jointly 
exploit the TOA information (between each target node and each reference node) and the RSS information (between 
neighboring target nodes) to accurately localize all the target nodes at once. The proposed strategy, thereafter referred 
to as cooperative TOA and RSS (COTAR), is advantageous in several ways: 

• COTAR operates on simple hardware and limited resources such as simple omni-directional antennas (rather 
than expensive directional antennas), and a narrow bandwidth (e.g. 2M Hz as in Zigbee sensors) rather than 
ultra-wide band. 

• Rather than treat sensors as independent nodes, COTAR makes clever use of the vicinity of two, three or more 
target nodes, and introduces among them very simple, practical and effective cooperation. The cooperation 
increases the signal spatial diversity and the geometric diversity, resulting in a significant improvement in the 
localization accuracy at the cost of only marginal overhead. 

• There has been consideration of jointly utilizing TOA and RSS in the literature [9][10], but the strategies thereof 
employ TOA and RSS in a rather natural and straight-forward manner, namely, each target node individually 
sends a beacon signal such that the M reference nodes record both RSS and TOA information and use all of 
the 2M measures to localize this target. By combining RSS with TOA, such hybrid schemes help improve 
the localization accuracy in short ranges, but as soon as the measurement range increases beyond 10 meters, 
they perform no better than TOA-only systems (yet with a higher complexity) [9] [10]. In comparison, COTAR 
combines TOA and RSS in a more clever manner and exploits each in its most favorable way - TOA for long 

'it should be noted that this type of error is different from the usual measurement error which can be modeled by a Gaussian random 
variable and which can in theory cancel out, given enough repeated measures. 
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ranges and RSS for short ranges, and hence achieves a significant performance gain in both short and long 
ranges. Additionally, COTAR can also effectively combat various channel conditions (e.g. fading, shadowing 
and NLOS) and measurement distortions, and is particularly suitable for heavily obstructed environments such 
as in an office building or in a forest. 

• The proposed strategy is rather general, allowing the joint localization of N (where N can be any number 
> 2) cooperative target nodes by exploiting M (where M can be any number > 3) reference nodes and two 
radio signal measurements (TOA and RSS). While the exact formulation takes a rather sophisticated, albeit 
systematic, nonlinear form, we have introduced a way to simplify it by approximating it to a low-complexity 
linear problem that can be iteratively refined in an efficient manner. Linearity makes incremental position 
tracking effective, and iterative refinement makes estimation accurate. 

• The proposed strategy is also analytically tractable. Asymptotic lower bounds on the localization distortion, 
quantified by the root mean square (RMS) distortion and the Cramer-Rao bound (CRB), are derived for COTAR, 
and simulations yield results extremely close to the theoretical bounds. 

Extensive simulations are conducted in static and mobile settings and in clear-sight and obstructed environments, 
and various practical issues and system parameters are evaluated, including the impact of the missing RSS informa- 
tion, the number of collaborating target nodes, the number and geometry of the reference nodes, the measurement 
range, and the mobility speed. Comparison to the conventional RSS, TOA, and TOA/RSS hybrid schemes confirm 
the superiority of the proposed COTAR scheme. 

The remainder of the paper is organized as follows. Section HIl introduces the system model. Section HID discusses 
the proposed COTAR localization protocol. Section |IV] formulates the maximum Ukelihood detection method and 
discusses its low-complexity linear approximation with iterative refinement. Section |V] analyzes the performance 
of the proposed strategy and presents a lower bound on root mean square error. Section |Vl] presents and discusses 
the simulation results. Finally Section IVIII concludes the paper. 



II. System Model 

Without loss of generality, we consider employing M (M > 3) reference nodes with pre-known coordinates vector 
(xr,yr) to help localize a cluster of target nodes with unknown coordinate vector (x, y) in a two-dimensional 
space. In the sequel, unless otherwise stated, we will use bold fonts (e.g. x) to denote vectors or matrices, and use 
regular fonts (e.g. x) to denote scalars. All the vectors are by default column vectors. Further, subscript s and r 
denote the quantities associated with the targets nodes and the reference nodes, respectively. 



A. Path-Loss Model and RSS Estimate Accuracy 

The average path loss of the radio signal propagation in general follows an exponential attenuation model, and 
can be expressed in unit of dB as: 

5(d) = 107?logio(d) +50, (1) 

where r] is the attenuation factor, d is the distance between the transmitter and the receiver, and go is the calibration 
pass-loss. The shadowing, multi-path and scattering phenomenons are collectively modeled as Gaussian distributed 
noise in dB (sometimes referred to as logarithmic Gaussian noise in the literature) with variance cr^. Thus the 
individual path-loss or signal attenuation follows a Gaussian distribution (when measured in dB): 

5(d)=N(5(d),a2). (2) 

Unless otherwise specified, we will use the attenuation factor -q = 3.086 and the standard deviation of the path-loss 
(Tg = 8 dB in the analysis and the simulations. 

From ([T]) and we can calculate the RSS estimation accuracy, measured by the standard derivation of the 
distance estimate based on RSS, as follows: 

std(d) = (3) 

This clearly indicates that the locahztion accuracy of RSS-based techniques increases linearly with the measurement 
range d. 
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B. Multi-Path Model and TOA Estimate Accuracy 

To study the TOA measurement error, consider modeling the multi-path channel profile as: 

h{T)=Y,ak5{T-Tk), (4) 

fc=0 

where Lp is the number of multi-path, and and au are the time delay and the attenuation coefficient of the /cth 
multi-path, respectively. To be practical, we assume that the receivers use a simple early-late structure and determine 
the signal arrival time by the first zero-crossing point of the correlator's output [14]. One major challenge about 
measuring TOA, especially in a highly obstructed environment with severe reflection and multi-path, is that the 
direct path (i.e. the first arrival path) is not necessarily the strongest path. As such, the receiver will likely miss 
the direct path, catch the strongest path and mistaken it for the direct path, thus arriving at a longer TOA estimate 
(and a longer distance estimate) than it actually i^ To account for this factor in our study, we use the two- 
ray approximation model developed in [14] to capture and simulate the effect of the strongest path. The two-ray 
approximation model consider two dominant paths: the direct path that occurs at some refrence time 0, and a second 
path that is usually a ground-reflected path (may be the strongest path) occurring after a small delay at time f . 
The small time dealy r is known as the mean excess delay [14]. The channel coefficient h is thus mathematically 
expressed as follows: 



Mr) = ^^<^(0) + ^^a<^(f), (5) 

where a is the Rayleigh faded random coefficient, and K is termed the Rician factor in the literature, which denotes 
the direct-path to multi-path energy ratio. A large K indicates that the channel is clear sight having strong LOS, 
while a small K indicates heavy obstruction and severe multi-path. 

Following Kim et. at.'s measurements in realistic engineering buildings [13], we consider two representative 
cases for the localization environment: the case of clear line-of-sight, which is associated with a mean excess delay 
f = 25.8 ns and a Rician factor K = 5, and the case of heavy obstruction, which has r = 76.9ns and K = 2. 
Additionally, following the performance predicting techniques in [14], we arrive at the the TOA estimation errors 
for different situations, which are listed in TableU and which will be used in our analysis and simulations. 

TABLE I 

T QA RSTTMATTON PRROR STANHARn nRBTVATTON FOR mFFPRFNlT T nr AT ,T7ATTnN FNVTRnNMFN TS 



situation 


Rician factor 


mean excess delay (ns) 


time-of-arrival error std (ns) 


clear line-of-sight 


5 


25.8 


8.8 


heavily obstructed 


2 


76.9 


40.2 



III. Cooperative TOA and RSS Localization Protocol 

The proposed COTAR strategy engages the set of target nodes in a simple but effective cooperation to facilitate 
each other's positioning. Consider the target nodes sending beacon packets one by one in a sequential manner. Due 
to the wireless broadcast advantage, when one target node transmits beacon signals, all the other target nodes in the 
vicinity can hear and measure the received-signal-strength from the transmitting node rather accurately. These RSS 
measurements contain valuable information about the relative distances and geological topology of the set of target 
nodes, which, when properly harnessed, can serve as useful calibration to compensate and mitigate localization 
errors. 

To exploit all this neighboring RSS information, we propose the following COTAR localization strategy depicted 
in Figure [T] 

The first target node starts out sending a beacon signal, and all the reference nodes measure and record the TOA 
of the beacon packet. At the same time, all the subsequent target nodes also measure the RSS from the transmitting 
target node. 

Starting from the second target node, each target node will embed all the RSS it has measured from the previous 
target nodes in the payload of its beacon packet, and broadcasts. That is, the second target node transmits to all 
the reference nodes the RSS information between the first target node and itself; the third target node transmits the 

^Note that 1 ns of TOA error amounts to 0.3 meter of distance error. 
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Reference 1 Reference 2 

Fig. 1. System model for COTAR with cooperating mobile noci^-knd M = 4 reference nodes. (A) Step 1: Target node 1 broadcasts a 
beacon packet, all the other (N — l) target nodes measure the RSS (from node 1), and all the M reference nodes measure the TOA (from 
node 1). (B) Step 2: Target node 2 embeds the RSS information it measured previously (i.e. RSS between node 1 and node 2) in its beacon 
packet and broadcasts, all the subsequent {N — 2) target nodes measure the RSS (from node 2), and all the AI reference nodes measure the 
TOA (from node 2). In addition, the AI reference nodes also try to decode the beacon packet to deduce the embedded RSS information. 
Next, target node 3 will air an RSS-bearing beacon packet. The procedure continues, until all the A'^ collaborating target nodes have sent a 
beacon packet. All the M reference nodes will combine all the TOA and RSS information and perform a joint position estimation for all 
the A'^ target nodes. 
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RSS information between the first target node and itself and between the second target node and itself; and so on. 
The M reference nodes not only measure and record the TOA of these beacon packets, but also try to demodulate 
and decode the beacon packets to obtain the embedded RSS information. 

Finally, the M reference nodes share all the TOA information (between targets nodes and reference nodes) and 
all the RSS information (between target nodes) they have collected and perform centralized computation, either at 
one of the reference nodes or a dedicated control center. The centralized computation estimates the positions of all 
the target nodes in a single batch, and may implement different estimation methods, some of which are optimal 
and others suboptimal. In this paper, we consider formulating the problem as as a maximum likelihood position 
estimation problem, and further propose a linear simpUfication of the ML method to render a fast and efficient 
implementation through iterative refinement. 

Here are a few comments: 

(1) The proposed cooperative strategy is simple and practical, and incurs very minimal cooperative overhead. 
Compared to non-cooperative localization strategies where the target nodes must each send a beacon packet anyway, 
the additional requirement imposed by COTAR is for the target nodes to record and include in its beacon packet the 
RSS measurements. In practical systems, any (wireless) packets, including load-less beacon packets, must consist a 
sufficient packet-head to provide such information as where the packet is originated from. In many systems, there 
is also a minimal payload length required for a wireless packet, such that the packet does not become too short to 
easily slip off without being detected. As such, in which case the RSS information actually gets a free ride. 

(2) The proposed strategy works for any set of M target nodes that are in the vicinity of each other and willing to 
cooperate, but a practical question arises as what is a good value for N. It is expected that the overall localization 
performance will improve with A'^ (with a diminishing return), but to engage many nodes in cooperation can be 
expensive or difficult in practical scenarios. Further, the communication overhead (recall that the ith target node is 
expected to convey (i — 1) RSS measurements to the reference nodes) as well as the complexity of the centralized 
estimation also increase with M. As we will show later in our analytical and simulation results, it is actually 
quite unnecessary to pull many nodes. Two to four cooperating nodes are sufficient, and in many cases, a mere 
engagement of two nodes (A^ = 2) promises to delivery most of the cooperative gain. 

(3) In the COTAR strategy, the reference nodes are expected to not only measure the TOA of beacon packets, but 
also extract from them the embedded RSS information. It is not necessary for all the reference nodes to correctly 
decode all the RSS information; as long as one out of the M reference nodes gets it (and conveys it to the central 
node), an optimal joint estimation can proceed. In the rare case when none of the reference nodes succeeds in 
deducing a particular RSS information, the system may adopt one of the following two strategies, (i) First, the 
M reference nodes can forward their individual reception of the beacon packets using, for example, amplify-and- 
forward (AF) or other ostensible signal-relaying technologies [15], to the central node for joint decoding. The 
central node is thus supplied with M noisy copies of the beacon packet with a diversity order of at least M, and 
has a considerably higher probability of extracting the RSS information than individual reference nodes can. (ii) 
Alternatively or as a last resort, the system can simply give up on the missing RSS information, and perform joint 
localization without having all the RSS information in place. Depending on how many and what RSS measurements 
are missing, the localization accuracy may be affected. We will discuss this issue in detail later, but in general, 
considering the rather small probabiUty that none of the M reference nodes has deduced the RSS information, the 
overall impact of missing RSS is very small. 

IV. Maximum Likelihood Position Estimation Using Joint TOA and RSS 

Having discussed the general COTAR protocol, we now elaborate the maximum UkeUhood algorithm the central 
node performs to simultaneously estimate all the target positions through the collection of TOA and RSS information. 

A. Path Loss among Target Nodes 

Let (x, y) be the horizontal and vertical coordinates of the actual positions for the targets nodes (x and y 
are each N dimensional column vectors). From the radio propagation attenuation model, we can write the average 
path loss between a pair of target nodes p and q, in the unit of dB, as: 




(6) 
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where go is the path-loss at the calibration distance, which is a constant and does not affect the locaUzation accuracy. 

The effect of multi-path, scattering, and shadowing phenomena may be collectively modeled as a multiplicative 
channel coefficient that follows a logarithmically Gaussian distribution. When expressing the channel attenuation 
in dB, the multiplicative coefficient becomes an addictive (noise) term Zi that follows the Gaussian distribution 
(in dB). Since there are target nodes, the actual path loss between any pair of target nodes can be expressed in 
an (^) dimensional column vector as 

B. Time-Of-Arrival from Target Nodes to Reference Nodes 

The distances between the ith reference node and all the N target nodes can be expressed using an M dimensional 
column vector function: 

di = ^(x,-x«)2 + (y,-y«)2, i = l...N. (8) 

The time of arrival (TOA) for the ith target node's signal arrives at all the M reference nodes will ideally take 
the form as follows: 

ti = di/c, ie[i..N], (9) 

where c is the speed of the waveform propagation, which is 299792458 m/s for radio frequency signals. We arrange 
all TOAs from each target nodes to each reference nodes into one single MN dimensional TOA vector t as 

t = [ti,t2,--- ,tjv]^. (10) 

All this TOA estimation error is modeled as an MN dimensional Gaussian-distributed error vector r. Hence, 
the noisy TOA vector from each target node to each reference node can be written as: 

t*=t + r. (11) 



C. Path Loss between Target Nodes and Reference nodes 
The average path loss between each target node i and all reference nodes, in the unit of dB, as: 

hi = lOr/logio (\J (x, - a;«)2 - (y, - y«)2^ +^3, g [l..A^],and p ^ q, (12) 

where qq is the path-loss at the calibration distance, which is a constant and does not affect the locahzation accuracy. 
We also rearrange all the RSS between each target node and each reference node in a single MN dimensional 
column vector h as: 

h = [hi,h2,--- ,hivF. (13) 

The effect of multi-path, scattering, and shadowing phenomena may be collectively modeled as a multiplicative 
channel coefficient that follows an addictive Gaussian noise term Z2 in logarithmical domain using unit dB. Since 
there are N target nodes, the actual path loss between each target node and any reference node can be expressed 
in an M X A?^ dimensional column vector as 

'^\mN)x1 — h(MAf)xl + Z2(MAr)xl- (14) 
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D. Unified RSS and TOA observation vector fiinction 

For a compact representation, we subsume the expected path-loss vector g and all the expected TOAs tj in a 
single column vector consisting of (^) + 2MN rows, termed observation vector: 

(") MN MN 

f(x,y) = (^,"t^,"h^f. (15) 
The true observation can be written as the sum of the observation vector f (x, y) and a Gaussian noise vector n: 

r=f(x,y) + n, (16) 

where the noise vector n is given by 

n = (Zi,r,Z2)^. (17) 

The observation error n denotes the measurement distortion caused by various imperfect conditions regarding RSS 
and TOA, and is assumed to be a multivariate random vector with D = (^) + 2MN dimensional positive-definite 
covariance matrix: 

A = E[{n- E[n]){n - E[n]f] 

( 2 ) MN MN 

= diag((7|^, (18) 



E. Simplification through Taylor Expansion 

The vector function g(x, y), t(x, y), and h(x, y) are nonlinear, which makes the measurement vector function 
f(x, y) technically challenging to handle. To make real-time computation possible, we approximate them using a 
single linear function by Taylor expanding them at the vector point ( xq, yo ). Mathematically, it can be written as 
follows: 

f(x,y)«f(xo,yo) + Gx (""'"'''Y , (19) 

where G is the partial derivative matrix for each measurement on each unknown variable. The partial derivative 
matrix can be further written a D x 2N dimensional matrix as follows: 



= [Gi,j]Dx2N 



dxi 
dxi 



dfn 
dxi 



X = XQ 

y=yo 



dfi 
dxj~i 

dh 

dXN 



dfn 

dXN 



X = XQ 

y=yo 



dyi 

dh 
dyi 



X=XQ 

y=yo 



dfn \ 
dyi I 



dh 
dyN 

dh 
dyN 



dfn I 
dyN I 



>C=3C0 

/=yo 



="0 



(20) 



whose elements jv] and Gj ,,g[jv+i,2Ar] represent the partial derivatives of the ith element in the observation 

vector with respect to the jth target node's horizontal coordinate .x*^^^ and vertical coordinate y'^^^ respectively. 
D = (^) + 2MN is the dimension of the observation vector f, and (xo,yo) is the initial position vector. 



F. Partial Derivative Matrix G 

A key computation step in the proposed joint estimation is the computation of the partial derivative matrix G, 
which consists of two parts: 




where the sub-matrices A is the partial derivative matrix pertaining to the RSS measurements between any two target 
nodes, B is the TOA measurements between target nodes and reference nodes), and C is the RSS measurement 
between each target node and each reference node, respectively. 
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The sub-matrix A has (^) rows and 2N columns, where rows correspond to pairs of target nodes, and columns 



Ai,2 
Ai,3 



correspond to the target nodes' horizontal and vertical positions. The exact formulation of A is defined as follows: 

/ Ai,2 \ 

(P / q) (22) 

A row vector Ap^g corresponds to the RSS between the pth and gth target nodes. It has four non-zero elements 
in the positions p, q, (N+p) and {N+q), and is zero everywhere else: 



/ dg 
' •••0, ^ 



,0---0 



V 




0---0, 



dg 



,0---0,— , , 

Qy{p) =c=xo Q {q) x=: 

•y v=vn ^ V = 





y=yn 



(23) 



pth qth {N+p)\h 

Plugging Q in (l23l) and after simplification, we obtain the following: 



{N+q)th 



dg 



dg 



dg 



X=XQ 

y=yo 



dx(i) 



Wt] ln(10)(x^^ 



X, 







y=yo 



(x, 



ip) 



X, 



dy(p) 

^ y = 



■0 

y=yo 



(24) 
(25) 



The partial derivative matrix for TOA, B, can be further decomposed to 2N sub-matrices, each corresponding 
to the TOA derivative matrix for one target nodes with respect to its horizontal and vertical coordinates: 



B 



/B.« 



V 



B 



(1) 



\ 



B 



(2) 



B 



(2) 



Ox 



B 



(N) 

y 



(26) 



Each submatrix B^,. is an M dimensional column matrix, denoting the partial derivatives of the M TOA values, 
measured by the M reference nodes about the zth target node, with respect to the ith target node's horizontal 
coordinates: 

T 



I 



b(^) 



9t 



(1) 



V 



dt\ 



(2) 



dt 



(M) 



' dx^^ 



(27) 



Submatrices By are similar, but the partial derivative is with respect to the the vertical coordinate pf the ith 
target node: 

T 



dt 



(1) 



dt 



(2) 



' dy^) 



dt 



{M) 



(28) 



Combining ^ and (O, and after simplification, the partial derivative sub-matrices can be written as follows: 



B« = - 







-y'o'r 










-xfY+{yf^ 


-yi^'r 




x(/^^-xi'> 





i G [l..iV] 



(29) 



Mxl 



9 



b(') 

y 









-y'o'Y 






y'o' 








Hyi"' 


-y'o'r 






-ylP 





i £ [1..N] 



(30) 



where c is the speed of Ught. 

The partial derivative matrix for RSS between each target node and each reference node, C, can be further 
decomposed to 2N sub-matrices, each corresponding to the TOA derivative matrix for one target nodes with 
respect to its horizontal and vertical coordinates: 



V 



.(1) 



\ 



i(2) 



(31) 



(i) 

Each submatrix Cx is an M dimensional column matrix, denoting the partial derivatives of the M RSS values, 
measured by the M reference nodes about the ith target node, with respect to the ith target node's horizontal 
coordinates: 

T 



dh 



(1) 



V 



dh 



(2) 



(M) 



.«=.(■) 



(32) 



(i) 

Submatrices Cy are similar, but the partial derivative is with respect to the the vertical coordinate pf the ith 
target node: 

T 



dh 



(1) 



dyii) 



dh 



(2) 



dh 



(M) 



.(.)=.« 



\ 



(33) 



Combining ([8]) and and after simplification, the partial derivative sub-matrices can be written as follows: 



a 



c(^) 



a 







-y^P 




ixl'^ 






-y'o'y 




y?'- 


-y'o' 








Hy','' 


-y'^'Y 






-y^'> 





i £ [1..N] (34) 



J\/xl 



Mxl 



, where a is the constant of lOryln(lO), and rj = 3.086 is the channel attenuation factor. 



G. Joint ML Estimation with Linear Approximation 

Having computed the partial derivative matrix G, we now discuss position estimation through maximum likelihood 
method and its efficient implementation. 

Since the noise in the observation vector function r(x, y) is assumed Gaussian distributed, we can write the 
probability distribution function (pdf) at the vector point (xo,yo) as: 

where |A| denotes the determinant of A, the superscript —1 denotes the matrix inverse, and D = (^) + MN is 
the dimension of A. 
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The maximum likelihood estimator calculates the value ( x, y ) that maximizes ( [35] ). Equivalently, it minimizes 
the quadratic term: 

Q(x,y) = [r-f(x,y)f A-i[r-f(x,y)]. (36) 

We further define 



ri = r - f(xo,yo) + G 



xo 

yo 



Plugging ( [37] ) into ( [36l ). the quadratic term (5(x,y) can be further written as: 



Q(x,y) 





T 











ri 



(37) 



(38) 



Minimizing Q(x, y) is a reasonable criterion for determination of an estimator even when the measurement 
error cannot be assumed to be Gaussian distributed. In this case, the resulting estimator is called the least squares 
estimator and A^^ is regarded as a matrix of weighting coefficients. Therefore, we turn the position estimation 
problem to solving the equation of 9(5(x, y)/9x9y = 0. 

From its definition, A is a symmetric matrix, i.e., A = A^. Since (A~^) = (A^) ^, it follows that (A^^) = 
A~^, which shows that A^^ is a symmetric matrix. Therefore, we can obtain: 



aQ(x,y)/axay = 2G^A-iG 



2G^A-iri = 0. 



We assume that the matrix of G^A ^G is nonsingular. Thus the solution of (|39l ) can be expressed as: 



G'^'A-^G) ^G^A-^ri 



(39) 



(40) 



Plugging (|37] ) in (l40l ). and after simplification, we arrive at the following estimate for the target nodes' position 
vector: 



+ (G^A-iG)-'G^A-i(r-f(xo,yo)) 



(41) 



where the first term denotes the initial positions of these target nodes, and the second term provides an incremental 
adjustment. 



H. Iterative Estimation Refinement 

Note that (|4TI ) suffers from the inaccuracy caused by approximating a nonlinear function f(x,y) using the linear 
Taylor expansion in ([T9l ) especially when (xo,yo) are far away from the actually position (x,y). To mitigate the 
error and improve the estimation accuracy in (|4T] ). we propose to solve this problem through iterative refinement, 
in which the estimated position of the previous iteration (x,y) is used as the initial position (xo,yo) for the next 
refining iteration. To start, when we have no a priori knowledge, we may take the center point (or any point) of 
the scenario as the initial location for all the target nodes. For example, suppose M = 4 reference nodes are placed 
in the corner of a square area with L meters per side, the proposed iterative localization proceeds as follows: 

^x[l] ■ 



y[i] 

m 
m 



(L/2,--- ,L/2) 



T 

2Afxl ' 



x[fc- 
Y[k- 



+ (G^A- 



^G) 



"Va- 



(42) 



(r-f(x[/c-l],y[fc-l])), A; = 2,3,. 



where [k] is the iteration index. 

It is expected that a good choice for the initial position (xo,yo) helps expedite the process while a bad one 
may delay the convergence, but this appears to be a non-issue, since our algorithm converges very fast. Through 
extensive simulations, we show that it generally takes no more than two to three iterations to arrive at an accurate 
estimation, even when the initial position is set very far from the actual position. 
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In the case of mobile tracking, when the previous-time position is used as the initial position to estimate the 

next-time position, a single iteration suffices. The iterative process not only improves the localization accuracy, 
but also makes the strategy particularly suitable for mobile localization. When target nodes are moving and their 
continuous trace needs to be detected, it is very natural to use their previous locations as the starting point to 
estimate the next locations. Unless the nodes are moving at extremely fast speeds, the new locations will not be 
far from the old ones, and hence a single iteration suffices, allowing efficient position tracking. 

V. Analysis of Estimation Accuracy and performance comparisons 

To provide a theoretical support, we first evaluate the Cramer-Rao lower bound for the proposed cooperative 
TOA-RSS localization scheme. We next analyze the theoretical root mean square error of the proposed maximum 
likelihood position estimation algorithm. We compare the bounds of the proposed localization scheme with the 
existing locahzation schemes. 



A. Cramer Rao lower bounds (CRB) analysis 

The CRB of an unbiased estimator (x, y) is given by 

cov(x,y) > /(x,y)-\ 

where is the Fisher information matrix (FIM) given by [17]: 

I{^,y) = -EV^,y (Vx,y^(r|x,y,Xr,yr)) , 

Ixy lyy 

where ^(r|x, y,Xr,yr) is the logarithm of the joint conditional probability density function. We can further write 

'"5^^(r|a:,y,xr,yr 



(43) 



T — —E 



Ixy — E 



"■yy 



-E 



0^^(r|a;,y,xi.,yr) 

dxdy 
'"a^^(r|x,y,Xr,yr) 



(44) 



The CRB on the standard derivation of each target node's position estimation is: 



(yCRB = . /min E [{x - x^ + {y - y)^] = i/mintr [cov(x,; 

= Vtr [I{x,yy^] 

/ Ixx ~l~ lyy 

Y Ixxlyy ~ Ixy 



(45) 



The expression of ^(r|x,2/,Xr,yr) is different for each localization estimation scheme. 
Existing Localization Schemes 

The traditional RSS scheme [18][19] uses only the RSS from a sensor node to all the M reference nodes (M 
RSS measures per sensor) to perform position estimation. We have: 



M 



(46) 



where ~ N (5(i(dB), al). 

Similarly, for the traditional TOA scheme [18] [19], only the TOA at the reference nodes M TOA measures per 
sensor is used for estimation. The logarithm of the joint conditional probability density function can be written as: 

M 

^TOA{r\x, y, Xr-, yr) = J] logPt,|x,y,a;.„j/,,, (47) 
i=l 
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where Pu\x,y,x„,y„'^^ {U, <7?), and Pu\x,y,x„,y„'^^ {U, o^r)- 

The existing non-cooperative hybrid TOA-RSS schemes exploit both the RSS and TOA from each senor to all 
the M reference nodes (, which exploits total 2M modality measures per sensor, including M modality measures 
per sensor for TOA and another M modality measures per sensor for TOA). As discussed in [17], for the hybrid 
TOA-RSS scheme, it is given by 

CTOA-RSs{r\x,y, Xr,yr) 
M M 

Pgi\x,y,Xri,y„ + ^OSPt,\x,y,x„,y,.i- 

(48) 

1=1 i=l 

We further substitute (|461 ) (|47] ) (|48] ) into (1431 ) (l44l ) respectively, and after some rather involved numerical calculation, 
we obtain the CRB for traditional RSS, TOA, and hybrid TOA-RSS localization schemes. 
The Proposed New Scheme 

The new cooperative TOA-RSS scheme developed here makes use of not only the TOA from each sensor to 
all the M reference nodes but also the RSS between the two neighboring sensors, to jointly estimate the poses of 
both sensors ( 2MN+{^) modality measures for N sensors). We can write the logarithm of the joint conditional 
probability density function as follows: 

M M N 

(49) 

1=1 i=l q=l,q^p 

where logPg^ jx,y ~ ^ {9p,q,(^l)- 

Similar to the procedure of calculation traditional localization schemes, substituting (l49l) in (l43l) and (|44l) and 
after simplification, we can compute the Cramer-Rao bound for the proposed new scheme. 



Traditional RSS scheme -c- Traditional TOA scheme 




Fig. 2. The CRB for different localization schemes at different locations in clear LOS situations. Four reference nodes are deployed on the 
four corners of a 18 by 18 meters' square. For the proposed COTAR scheme, 4 target nodes cooperates localizing for each other in a 1 by 
1 meter's grid. 

The Cramer-Rao bounds for different schemes are shown in Fig|2l where 4 reference nodes are deployed at the 
four corners of an 18 x 18 meters' square scenario in clear line-of-sight situation. Four cooperative target nodes 
help each other in the proposed COTAR scheme at the four comers of a 1 x 1 meters' square. It is clear to see that 
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the RSS only scheme has a poor localization accuracy from 8 to 10 meters. The traditional TOA scheme obtains 
a better accuracy of 2.7 meters in the middle and 2.9 meters at the four corners. The hybrid TOA-RSS scheme 
helps improve the accuracy of traditional TOA scheme at the four corners to 2.2 meters, but not in the middle of 
the scenario. The proposed scheme has an excellent localization accuracy, which achieves 1.55 meters' accuracy 
in the middle and 1.2 meters' accuracy at the four comers of the scenario. 



RSS 

TOA (clear LOS) 
hybrid (clear LOS) 
new (clear LOS) 
TOA (heavily obstructed) 
hybrid (heavily obstructed) 
new (heavily obstructed) 




lOOOOOOOOOOOOOOOOOOOOOOOCK > 



100 

scenario size (meter) 



200 



Fig. 3. The average CRB for different localization schemes vs increase of scenario size. Four reference nodes are deployed on the four 
comers of a increasing size's square. For the proposed COTAR scheme, four target nodes cooperate localizing for each other in a 1 by 1 
meter's grid. Both clear line-of-sight and heavily obstructed situation are considered. 

We further invest the localization accuracy vs increase of scenarios size for different localization schemes in 
different channel conditions in Figl3l It is shown that the RSS scheme's localization error increases linearly fast, 
which makes it useless in a large scenario. On the contrary, the accuracy of traditional TOA scheme does change 
with the increase of range. The hybrid TOA-RSS scheme can only improve the accuracy when the scenario size 
is less than 10 meters' in clear line-of-sight situation, and 100 meters in heavily obstructed situation. Otherwise, it 
can only obtains the same performance as traditional TOA scheme. For the proposed COTAR scheme, it performs 
excellently in any scenario size. It cuts half down the localization error from 12 to 6.5 meters in heavily obstructed 
scenario, and improves the accuracy from 2.7 meters to 1.6 meters in clear line-of-sight situation. 



B. Root Mean Square Error Analysis for ML Estimation 

We now analyze the performance error of the proposed detector in (|4TI ). and evaluate the errors with respect to 
different positions (x, y). 

The expression in (|4TI ) shows that when the measurement vector r is Gaussian distributed, the estimation vector 
becomes a 2N dimensional Gaussian random vector with a pdf: 



T 



p-1 f^-E 



(2^)^|P|i/2 



(50) 
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where P, the covariance matrix of the estimation vector ( ^ ), is defined as: 



E 



E 



E 



(51) 



There are two types of inaccuracy associated with a Gaussian distributed variable, the mean bias p and the variance 
ip. The estimation error e should in general consist of both: 



p + tp. 



S) (y 

On the other hand, substituting ([T6l) into (|4TI ) and rearranging terms, we can write the estimation vector as: 



(52) 



X 



X 



f(x,y)-f(a,b) 



X-Xo' 


+ n 




.y-yo. 





(53) 



It is clear from ( 153] ) that the estimation error comprises both a linearizing error and a noise term. The linearizing 
error in (1531) contributes to the mean bias: 



P-- 



G^A-iG)-iG^A-M f (x, y) - f (xq, yo) - G 



X-Xo 

.y-yo. 



(54) 



The estimator will be unbiased, if f (x, y) is a linear vector function; but since it is not, the inaccuracy of the Taylor 
expansion in ( fT9l ) thus introduces non-zero mean bias to the estimation. From ( fT9l ). we also know that if the initial 
position (xo,yo) is sufficiently close to the actual position xq — )■ x and yo — )■ y, the linearizing error p in (|54l ) 
will become vanishingly small and can be safely ignored. 

The residual noise in the estimation, caused by the measurement error vector n, results in variation in the 
estimation results. This noise part is expressed as: 

fG^A-iG)-iG^A-i' 



n 



(55) 



=xo,y=yo- 

Here the error vector n subsumes all the error contributors, including the measurement distortion of the signal TOA 
and the RSS between the two cooperative nodes, and others uncertainties in the system. It is reasonable to assume 
that the measurement distortion n for TOA and RSS are unbiased, such that ip is, & 2N dimensional column vector 
following a zero-mean (vector) Gaussian distribution. 

Since the mean bias is a constant vector given the initial position 
P of the estimation error is only affected by ( [55] ): 



yo 



the calculation for the covariance matrix 



E 



(G^A 



-E[V])(V- 



my 



(56) 



Note that the diagonal elements of P denote the variances of the errors in the estimated positions. 

Since the measurement distortion n is Gaussian distributed with zero mean, the maximum likelihood or the least 
square estimator for the linearized model is the same as the minimum variance unbiased estimator. A scalar measure 
of the estimator accuracy is the root mean square error e, which is defined as 



v/tr(P)/iV. 



(57) 



The RMS error bound of the proposed ML estimator for the COTAR scheme is plotted in Fi^ as a function of 
the distance between target nodes and the number of cooperative target nodes. For comparison, we also shown the 
CRB of the COTAR scheme in the same picture. We consider four reference nodes are deployed at the four corners 
of a 1000 X 1000 meters' square scenario, and target nodes are placed in a (\/]V A) x (\/iV A) meters' grid. The grid 
size A is the distance between two closest cooperative target nodes. It is well-known that the CRB sets the lower 
bound of a scheme for all the possible detectors, and may be achieved by an unbiased optimal detector. That our 
root mean square error bounds matches perfectly with the CRB (shown in Fig. HJ), and that our simplified iterative 
algorithm practically achieves these bounds (as will be shown in simulations), clearly indicates the efficiency of 
the proposed detection algorithm. We also find that the gird size A and the number cooperative target nodes are 
two important factors to the localization accuracy. For example, increasing the number of cooperating nodes from 
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Fig. 4. Performance comparison between joint maximum likelihood estimation algorittim vs Cramer Rao lower bound for the proposed 
COTAR scheme under different number of cooperating nodes and distance. 

1 (non-cooperative) to 25 can effectively improve the localization accuracy from 2.7 meters to 0.7 meters. It also 
shows that the cooperative target nodes should not be too far away. Due to the quickly-decreasing accuracy of RSS 
with distance, a close neighbor can do more to improve localization accuracy than several remote collaborators. 



A. COTAR vs Existing Strategies 

We now conduct computer simulations to study the realistic performance for the proposed COTAR localization 
algorithm. We firstly consider four reference nodes (M = 4) are located at the four comers of a square area with 
edge length L = 50 meters. 

To start, we consider four cooperating target nodes (A^ = 4) that move together in the square area and always 
stay at the four corners of a 1 x 1 meter's grid. Four localization strategies based on RSS, TOA, hybrid TOA-RSS, 
and the proposed COTAR, are evaluated: 

1) RSS-only: Each target node sends a beacon packet independently, and all the M reference nodes measure the 
RSS information, and route them to a central node for maximum likelihood estimation. The target nodes do 
not cooperate, and M RSS measures are used to estimate the position of any one target node. 

2) TOA-only: It is similar to RSS-only, but the M reference nodes measure the TOA information instead. M 
TOA measures are used to estimate the position of any one target node. 

3) Hybrid TOA/RSSE [9][10]: It is similar to RSS-only and TOA-only, but the M reference nodes measure both 
the TOA and the RSS information. The scheme is again non-cooperative, and 2M measures (M TOA and 
M RSS) are used to estimate the position of any one target node. 

4) COTAR: As discussed before, the scheme is one that combines TOA and neighboring-RSS. A set of N 
target nodes are engaged in a simple cooperation, and MN TOA measures, (^) neighboring-RSS measure, 
and NM remote RSS information are used to jointly estimate the positions of the target nodes. In the 

^For convenience, we call the RSS between target nodes and reference nodes "remote-RSS" and the RSS between two (nearby) target 
nodes "neighboring-RSS." 



VI. Simulation results 
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Fig. 5. Impact of the number of cooperating nodes on the performance of COTAR for different Rician distribution factor K. 50 x 50 
meters grid, M = 4 reference nodes at the comers, bandwidth=2Mhz. 



simulation, the scenario center with coordinates is always used as the initial position for both target nodes, 

and two iterations are used in the iterative refinement. 
The localization results for tracking mobile nodes are shown in Fig|5l where red zigzag line shows the moving 
track of a group of four cooperative target nodes and blue dots represents the estimated positions of different 
schemes. The results of 15 experiments are showrQ and in each experiment, position tracking is performed in a 
continuous manner, with the previous estimated location serving as the starting point for the next one, and two 
iterations to estimate the incremental errors. Also shown in the figure are the performance contour. Because of the 
large distance between each target nodes and the reference nodes, we are not surprised to see that the RSS-only 
method performs rather poorly, with the estimated positions scatter all over the place, yielding an estimation error 
RMS of as large as 10 to 12 meters. In comparison, the traditional TOA-only method is much less sensitive to the 
distance. It attains a satisfactory performance of 2.7 meters of error. When TOA and remote-RSS are combined as 
hybrid TOA-RSS scheme [9], only the performance of the four corners of the scenario is effectively improved(, 
where reference nodes are deployed), but not the most of other positions especially in the middle of the scenario. 
The proposed COTAR strategy provides a much better solution by leveraging both the good accuracy of the RSS 
between two neighboring target nodes in short distance and the large localization range from TOA techniques. 
COTAR achieves a high localization accuracy of less than 1.5 meter's of distortion in the middle of the scenario, 
and 1.2 meters' accuracy at the four corners. 

B. Impact of Initial Position and Number of Iterations 

We next study the impact of the number of iterations and the initial position on the localization accuracy. The 
same 50 x 50 meters' square area with M = 4 reference nodes and = 4 cooperating target nodes is simulated. 

''We performed more experiments, but they make the picture too dense to see. 
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We consider the clear line-of-sight channel situation, the cooperative grid size to be 1 meter, and use the scenario 
center as the initial position for the target nodes. 

We evaluate the gap between the actual localization error and the theoretical lower bound. This gap is plotted 
in Figure OA) and Figure OB) for one and two estimation iterations, respectively. We see that the accuracy of 
the starting point can make a difference, but only to the first iteration. As shown in Figure OA), with only one 
iteration, the localization strategy may achieve a tiny RMS error of < 0.01 meter from the lower bound near the 
grid center, but an accuracy gap of as large as 1.9 meters near the four edges. It is interesting to observe that, 
when the target nodes get very close to the corner points, their closeness to one of the reference nodes, and hence 
the accurate TOA measures, can help mitigate the negative impact caused by the inaccurate starting point (and the 
insufficient iterations). As the iterations increases to two, the accuracy gap quickly drops to near-zero across the 
entire region. In the simulated 50 x 50 (meter) grid, it appears that when the initial position is not offset by more 
than 10 meters from the actual target locations, one iteration is enough; Otherwise, two iterations suffice to bring 
down the estimation error to the lower bound. 

C. Impact of Number of Reference Nodes and Cooperative Nodes 

We first consider there are more than 4 reference nodes, say 9 reference nodes in the scenario. The 9 reference 
nodes are evenly deploy at every 25 meters as black ■ in FigjT] We study the performance of the hybrid TOA-RSS 
scheme and the proposed COTAR scheme and plot the simulation results in Fig|7] From this figure, we see that the 
ML estimation algorithm works robustly in any position in the scenario, and increasing reference nodes number 
from 4 to 9 can effectively improve the localization accuracy for both hybrid TOA/RSS and COTAR schemes. At 
the middle of the scenario, the localization accuracy improve from 2.65 to 1.65 meters for hybrid TOA-RSS, and 
1.7 to 1.05 meters for COTAR. 

We further study the localization performance with more than 4 cooperative nodes. The same system set up is 
used, and we consider the cases when N = 9andl6 cooperative target nodes are available for cooperation. 

The simulation results, plotted in Figure [8j clearly speaks for the importance of increasing number of cooperative 
nodes. With the help of 9 cooperating target nodes, COTAR can achieves 1.05 meters' accuracy in most of the 
area. If we further increase the number of cooperative nodes to 16, we are glad to find that COTAR obtains an 
excellent cooperative localization performance of only 0.85 meters' error. 

D. Impact of Missing RSS 

As mentioned before, although uncommon, it is possible for some RSS information to miss in the equation, either 
because a target node fails to estimate it in the first place, or because all the M reference nodes fail to decode it 
from the beacon packet. In such cases, the joint estimation should proceed without this information. One possible 
way to handle such abnormality is to reformat all the relevant matrices in (l22l) by extracting the con^esponding 
row(s) from (l23l) . A more convenient way is to simply set the missing RSS information to be zero in the original 
formulation in ( [231) without changing the dimension of the partial derivative matrix. 

The potential impact of missing RSS information on the COTAR performance is presented in Figure |9l Cases with 
different cooperating nodes having different probability (percentage) of missing RSS are evaluated. It is encouraging 
to see that the localization accuracy degrades rather gracefully with missing RSS. When the missing RSS is below 
20 percentage, the increase in average RMS error is very minor - no more than 0.07 meters for 2 cooperating nodes 
and negligible for 16 cooperating nodes. When all the RSS is missing (100% missing probability), the COTAR 
strategy reduces to the non-cooperative TOA-only scheme with RMS error of about 2.67 meters. It is also evident 
from the simulation that increasing the number of cooperating nodes is an effective way to compensate RSS loss 
and reduce localization error. 

E. Tracking Highly Mobile Nodes 

The proposed localization strategy is also evaluated in a highly mobile system such as the 3G mobile cellular 
phone networks. We consider four base stations situated at the four comers of a square grid with edge length 
L = 1000 meters, acting as the reference nodes {M = 4). Suppose a pair of cell phones are moving together in 
the grid with a high speed of up to 160 km/hour. Suppose the cell phones move with random directions and when 
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Fig. 6. (A) RMS practice-theory gap (difference between the simulated RMS errors and the theoretical lower bound) with only one iteration. 
(B) RMS practice-theory gap with two iterations. 50 x 50 meters grid, M = 4 reference nodes at the corners, N = 2 cooperating target 
nodes, bandwidth=2Mhz, Rician factor K — 5. 
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Fig. 7. Impact of increasing the number of reference nodes to 9 for the performance of COTAR. 50 x 50 meters square scenario, M = 4 
reference nodes at the comers, bandwidth=2Mhz. 
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Fig. 8. Impact of increasing the number of cooperating nodes on the performance of COTAR, 50 x 50 meters square scenario, M — 4 
reference nodes at the comers, clear line-of-sight. 



they hit the boundary, they bounce off like the reflection of a ray of light. The COTAR algorithm is used to track 
the mobile positions of these cell phones, where the previous estimate is used as the starting point for the next 
estimate, and either one or two iterations of refinement is performed. To make the system simple and light-load, 
we consider a very low sample rate, such that beacon packets are transmitted and localization is performed once 
every 5 seconds. 

The localization RMS error is plotted at a function of the mobile speed in Figure [TO] We can see that when 
the nodes move at a slow to medium speed (< 80 km/hour), 1 iteration is enough to obtain good accuracy. As 
the mobile peed increases beyond 80 km/hour, since our sample rate is rather low and two succeeding sampling 
points are fairly separated, and since the localization errors may accumulate, 1 iteration becomes a little insufficient, 
resulting an increase of localization RMS error from 2.55 meters (at 80 km/hour) to 3.60 meters (at 160 km/hour). 
To bring down the error, one can simply add one more iteration, such that the localization error is kept consistently 
low and close to the lower bound. 

VII. Conclusion 

Low-cost and accurate wireless localization is challenging, due to the limited bandwidth, hardware and computa- 
tional resources, as well as possibly severe multi-path fading and non-LOS conditions. In the paper, we propose an 
effective cooperative-TOA-and-RSS localization strategy which is simple, practical, and inexpensive to deploy, and 
which requires no additional hardware than the usual cheap sensors (like Zigbee nodes). By engaging two or more 
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Fig. 9. Average localization accuracy as a function of the percentage of missing (neighboring) RSS information. 50 x 50 meters scenario, 
M = 4 reference nodes at the four corners, A'^ = 1, 2, 4, 9, 16, 25 cooperating nodes, clear line-of-sight, target grid size to be 1 meter. 




Fig. 10. Localization accuracy for mobile nodes at different moving speeds and iteration numbers. 1000 x 1000 meters grid, M — 4 
reference nodes at the corners, N — 2 cooperating mobile nodes which are 1 meters' apart, clear line-of-sight situation. 
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target nodes in simple but effective cooperation, and by leveraging useful estimation and detection techniques, the 
new COTAR strategy reaps the merits of both TOA (i.e. decent performance fora a large localization range) and RSS 
(i.e. high accuracy in short distances), without incurring much cooperation overhead to the target nodes. Analytical 
and simulation results confirm the good performance of the new strategy in both clear and heavily obstructed cases. 
The proposed COTAR strategy is particularly useful for (high-speed) mobile locaUzation, since the algorithm can 
naturally and efficiently achieve incremental update through iterative refinement. 
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